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ABSTRACT 

We consider (l+l)-dimensional QCD coupled to scalars in the adjoint representation 
of the gauge group SU(A^). This model results from dimensional reduction of the (2+1)- 
dimensional pure glue theory. In the large-iV limit we study the spectrum of glueballs 
numerically, using the discretized light-cone quantization. We find a discrete spectrum of 
bound states, with the density of levels growing approximately exponentially with the mass. 
A few low-lying states are very close to being eigenstates of the parton number, and their 
masses can be accurately calculated by truncated diagonalizations. 
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1. INTRODUCTION 



The light-cone quantization supplemented with a regulator in the form of discretized 
longitudinal momenta is a promising approach to nonperturbative QCD [1,2,3], which may 
provide tools for calculating the hadron spectrum as well as the interaction cross-sections 
starting from first principles. It has been successfully used in studies of various field theories 
in 1-1-1 dimensions, and some progress has been made in generalizing it to (3+l)-dimensional 
theories [2,4]. Recently, Dalley and one of the authors applied this method to a new type 
of models, namely, two-dimensional large- QCD coupled to matter in the adjoint repre- 
sentation of SU(iV) [5]. These models are far more complex than the large- QCD coupled 
to matter in the fundamental representation [6-9], but they are the simplest models where 
one can study some genuine QCD effects. There are two color flux tubes attached to each 
quantum and in that sense the quanta of the adjoint matter resemble physical gluons. The 
glueball-like bound states may contain any number of quanta connected into a closed string 
by the color flux tubes. Pair creation inside the string is not suppressed in the large- 
limit, leading to eigenstates which are mixtures of strings with different numbers of partons. 
(In the 't Hooft model, due to the absence of transverse gluons, one finds only the quark- 
antiquark bound states [6].) Therefore, the adjoint matter imitates some transverse gluon 
effects. 

Another way to see this is to dimensionally reduce (2-|-l)-dimensional pure glue theory, 
and to note that the zero mode of the transverse gluon field acts as the adjoint matter field 
coupled to (1-1-1) -dimensional QCD. Consider the pure glue SU(A^) theory in 2-1-1 dimensions. 



Compactifying one of the spatial dimensions, y ~ y + L, we may ignore the dependence of 
fields on y as L ^ 0, i.e. dAi^/dy = 0. The action then reduces to 



where — g1/L, and 0(a;°, x^) — Ay/g is a traceless N x N hermitian matrix field, whose 
covariant derivative is given by Da4> = da4> + ^[^a, 4']- Therefore, (p represents the remnants 
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of the transverse gluon degrees of freedom. In order to absorb the logarithmically divergent 
mass renormalization it is necessary to add a mass term for (p: 



The light-cone quantization and the spectrum of the theory (3) were considered in ref. [5]. 
In this paper we present results of a new numerical diagonalization. We find that holding 
the renormalized mass fixed results in a convergent and physically reasonable spectrum of 
glueballs. We compute this spectrum by solving the linear light-cone Schrodinger equation 
describing the bound states. The density of states grows roughly exponentially with their 
mass, as expected in a physically interesting gauge theory. We also find that a few low-lying 
states are very close to being eigenstates of the parton number, while a typical excited state 
is a complex mixture of states with different parton numbers. These results are similar to the 
results of our previous work [10] where we performed a detailed numerical study of a similar 
model, with the adjoint scalar replaced by an adjoint Majorana fermion. For a different 
approach to the model with the adjoint fermions see ref. [11]. 



In this section we summarize the light-cone quantization of (1 + l)-dimensional SU(A^) 
gauge theory coupled to a scalar field in the adjoint representation [5]. Consider the action 
(3) and choose the light-cone gauge A- = 0. Introducing the light-cone coordinates — 
(x^ ± x^)/y/2, we obtain 



where — i[0,9_0]y is the longitudinal momentum current. A similar model, with 
taken to be a general complex matrix, was considered in ref. [12] in an attempt to solve 
higher-dimensional pure glue theory. 




(3) 



2. LIGHT-CONE QUANTIZATION 
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In the light-cone quantization is treated as the time and the canonical commutation 
relations are imposed at equal a;+: 

[<Pij{x~), d-(f)ki{x~)] = ^ 5{x~ - x~)5ii5jk . (5) 

The action does not depend on the time derivatives of the gauge potential which, there- 
fore, can be eliminated by its constraint equations. The light-cone momentum and energy, 
— J dx~T~^^, are found to be [5]: 



P+ = y dx- Tr[(a_0)2] , 



(6) 



Our goal is to solve the eigenvalue problem 

2P+p-|$) = m2|$) , (7) 

where the physical states must satisfy the zero-charge constraint 

Jdx-J+\^)^0, (8) 

arising from integration over the zero- mode of Since [P~^, P~] = 0, 1$) is a simultaneous 
eigenstatc of P+ and P^. In practice it is easy to ensure that |$) carries a definite P"*", but 
the subsequent solution of eq. (7) is highly non-trivial. 

Introducing the mode expansion 

oo 

<Pij -^j^ (a,,-(fc+)e-^^"-- + a],{k^y^'^-) , (9) 



with 

.+^ (1+ 



,(A;+),a;fc(fc+)] = - , (10) 
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eq. (6) can be written in terms of the oscillators as follows: 



oo 

P+ = jdk ka\.{k)aij{k) , (11) 



p-^^rn'^J ^al{k)aijik) 





oo 



+ ^^^^^'^^^^ =^|^'^(^i + k2-k3- k^)a\.{ki)a}j^{kA)aki{ki)aii{k2) 



+ B5{ki + k2 + k3- k4){alj{ki)a'ji{k2)al^{k3)aki{k4) + alj{k4)aki{ki)aii{k2)aij{k3))^ , 

(12) 

where 

_ (^2 - ki){k4 - fcs) _ (^3 + ki){k4 + k2) 
{ki + k2Y {kA-k2Y ' 

(13) 

{ki + kA){k3 - k2) {ks + kA){ki - k2) 
{k3 + k2f ^ {ki + k2Y 

and we have dropped the superscript plus on ki for brevity. The renormalized mass m is 
defined by 

oo k 

2 2^9^N [dp g^N f k 



Note that the second term in the mass renormalization is finite in the principal value sense, 
while the first term is logarithmically divergent. We will see that the spectrum is finite if 

the renormalized mass m is held fixed. 

An important simplification in the light-cone quantization is that the oscillator vacuum 
satisfies 

P+|0) = 0; P-|0) = 0. (15) 

Other Fock states are then constructed by acting with creation operators ajj on the vacuum. 
The zero-charge condition (8) requires that all the color indices be contracted. Therefore, 
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we look for eigenstates of eq. (7) in the form 

P+ 

00 » J 

j=l «=i 

fjih, k2,..., kj)N-^/^ Tr [a^ki) . . . a^kj)]\0) , 
where the wave functions have cychc symmetry 

fi{k2, k3,..., ki, ki) = fi{ki, k2, ks, . . . , h) . (17) 

This state is trivially an eigenstate of , and the problem is to ensure that it is an eigenstate 
of P-. 

The complexity of the coupling to adjoint matter arises mainly from the presence of pair 
production and pair annihilation terms in P~ . These terms appear in the leading order 
of the expansion, provided that g'^N is kept fixed in the large- limit. On the other 
hand, the terms in P~ that take one color singlet into two are suppressed by 1/N , which 
implies the stabihty of bound states in the large- A?" hmit. The glueball wave functions satisfy 
linear eigenvalue equations [13, 5], which, upon introducing longitudinal momentum fractions 
Xi — kf / P'^ ^ can be written as 

X2, . . . , Xi) = ^ ^ fi{xi,X2, . . . , Xi) + 4^^= ^2, . . . , Xi) 

X1+X2 

+ J dy xi+X2-y,X2„..., Xi) {x2 - xi) {xi + X2 - 2y) 



^\/xiX2y{xi + X2-y) ixi + X2)'^ 

f dy {xi+y){xi + 2x2-y) ^ \ ff , m 

+ / 7 To , / , = [fi{xi,X2,...,Xi) - fi{y,xi + X2-y,X3,...,Xi)\ 

J [xi - yY A^xiX2y{xi +X2-y) 





X1+X2 



+ ld I ~ ^ ~ ^'^2, • • • ,a;i) 



^ ^ 4:y/xiyz{xi - y- z) 



(xi + y){xi -y-2z) ^ (2xi - y- z){y - z)' 



+ 



(xi - |/)2 {y + z) 

/i-2(xi +X2+ X3, X4, . . . , Xi) r(2xi X2 X3)(x3 - X2) , (xi X2 2x3)(xi - X2) 



+ 



4v/xiX2X3(xi +X2 + X3) L (X3+X2)2 (xi + X2)^ 

± cyclic permutations of (xi, X2, . . . , Xj) , 



(18) 
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where = 2P^P is the bound state (mass)^. In writing this equation we have used a 
principal value integral [8] 



dy 



{xi + y){xi + 2x2 - y) 



— — TT . 



(19) 



^y{xi + x2 - y) 







Note that the only singular term in eq. (18) is well-defined in the principal value sense. 

By truncating eq. (18) to the two-body sector [see eq. (23)], we can show using a constant 
variational wave function that is unbounded from below for < 0. For m? > 0, on 
the other hand, the of every state is positive. Thus, the critical value of w? is zero. 
We believe that this is the correct value to consider if we want to regard our model as the 
proper dimensional reduction of the (2-1-1) -dimensional pure glue QCD. Therefore, we will 
be interested primarily in the — case. 

The form of eq. (18) indicates that the interactions conserve the number of partons 

modulo 2. Therefore, there are two kinds of bound states: those involving mixtures of states 
with even parton numbers, and those involving mixtures of states with odd parton numbers. 
From here on, we will call them even and odd glueballs respectively. Furthermore, eq. (18) 
possesses another Z2 symmetry, T [11]: 



The Z2 quantum number T has two possible values, 1 and —1. In terms of the original 
field, T : (f)ij — > which obviously leaves invariant [11]. Physically, every bound state 
can be thought of as a superposition of oriented closed strings, and the quantum number T 
describes the transformation property under a reversal of orientation. 

3. THE DISCRETIZED APPROXIMATION AND NUMERICAL RESULTS 

The system of equations (18) involves an infinite number of multivariable functions. It 
seems very difficult, if not impossible, to solve this system of equations analytically. There- 
fore, here we follow the same strategy as in [10] and replace the continuum equations (18) 
by a sequence of discretized approximations, such that the eigenvalues of the discretized 



fi{xi, X2,..., Xi) = Tfi{xi, ...,X2, xi) . 



(20) 
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problems eventually converge to the eigenvalues of (18). In the light-cone quantization, a 
simple discretized approximation is obtained by replacing the continuous momentum frac- 
tions a; by a discrete set n/ where n are odd positive integers, and the positive integer K is 
sent to infinity as the cut-off is removed [3,2]. The restriction to odd integers corresponds to 
the choice of antiperiodic boundary conditions in the discretized light-cone quantization, and 
leads to the best convergence towards continuum limit. Thus, the functions fi{xi, X2, . . . ,Xi) 
are replaced by 

ffiK, n.) = if <'-"V.(|,|,...,^) 

with Y^l^i rij — K, and 



i 

/ 



2 ^ 

odd n>0 



For the even glueballs, gi{ni, n2, ■ ■ ■ ,ni) are normalized as 



i-l 



5Z ••• \9i{'ni:n2,...,ni-i,K -^rij)^ =1, (21) 

even i odd ni>0 odd ni_i>0 j=l 

and similarly for the odd glueballs. 

The number of independent components of gi{ni,n2, ■ ■ ■ ,ni) is equal to the number of 
partitions of K into positive odd integers, modulo cyclic permutations. Finding all such 
partitions is a combinatorial problem that is easily solved with a computer program. If we 
regard all the independent gi{ni,n2, . . . ,ni) as components of a vector, then the discretized 
eigenvalue problem (18) is equivalent to diagonalizing a matrix 

M^TT 

-^^KixV + T), (22) 

2 

where x — is the dimensionless parameter. The main difficulty in this approach is that 
the number of states increases rapidly with K. For example, there are 765, 1169, 1810, 
2786, 4340, 6712 states for K = 20,21,22,23,24,25, respectively. In actual calculations it 
is advantageous to consider separately the even and odd sectors under T, and introduce the 
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linear combinations 



gf{ni,n2, . . . , nj) = (^gi{ni, n2, . . . , Ui) ± giirii, Ui-i, ni) j . 

The entries of V and T can be read off from eqs. (18). For example, for K = 6 there are 4 
partitions which correspond to the 22-even sector: {(1, 1, 1, 1, 1, 1), (3, 1, 1, 1), (5, 1), (3, 3)}. 
The matrix to be diagonalized in this case reads: 



K{xV + T) 



/36a; + 97r ^ \ 

20x+(3+V3)7r 

n rv 36x I 16+37r / 512 

^ ^ 5 + V5 V 5 

V 4x+^ + 7ry 



For this simple example there are no partitions corresponding to the Z2-odd sector. For 
larger values of however, the total number of partitions splits approximately equally 
between the even and odd sectors, allowing us to reach higher values of K. For K — 2A 
there are 2358 components in the Z2-even sector and 1982 components in the Z2-odd sector. 
For K — 2b there are 3544 and 3168 even and odd components, respectively. 

A good numerical procedure is to calculate the spectrum for a fixed x and a range of 
values of K, and then to extrapolate the results to infinite the continuum limit. We will 
also assume that some bulk properties of the spectrum can be estimated from the results at a 
fixed large K. We will be most interested in x = which corresponds to the limit of massless 
quanta. Note that, with our choice of the regulator, we observe only even (odd) glueballs 
for an even (odd) K. We will present here the results for our biggest diagonalization in each 
sector separately: for even K — 2A and for odd K — 25. Our analysis is very similar to the 
one performed in ref. [10] for the model with the adjoint fermions. 

In fig. 1 we show the spectrum of even glueballs for x = and K = 24, with the mass 
plotted vs. the expectation value of the number of partons, n. The density of states obviously 
increases rapidly with the mass, and almost all the states lie within a band bounded by two 
(n) ~ M lines. One interesting feature of our results, already noted in ref. [5] for smaller 
K, is that for a few low-lying eigenstates the wave functions are strongly peaked on states 
with a definite number of partons. For example, for K — 2A the ground state has probability 
0.99899 to consist of two partons, and the first excited state has probability 0.99734 to consist 
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of four partons. For the ground state, it is a good approximation, therefore, to truncate eq. 
(18) to the two-body sector described by the wave function (j){x) = f2{x, 1 — x). From eq. 
(17) we have 4>{x) = 4>{1 — x). Thus, we obtain 

= + + ^ / dy^^^tl^ {x + y){2-x-y) 

^' ^(1-^) 2 ^x{l-x) 2n J {x-vY ^x{l - x)y{l - y) 

(23) 

This is simply the bound state equation for the theory with scalar quarks in the fundamental 
representation of SU(iV) [8,9], with — > 2g^. The interaction strength is doubled because 
now there are two flux tubes connecting the pair of partons. The lowest eigenvalue of eq. 
(23), M2 = 4.76c/27V/7r, provides a good upper bound on the lowest eigenvalue of eq. (18). 

As the excitation number increases, however, the wave functions typically become quan- 
tum superpositions of states with different parton numbers. The physical picture is that 
a typical excited state contains some number of virtual pairs, and our data supports this 
expectation. In order to quantify this effect, we will call a state pure if it has probability 
> 0.9 to be in one of the number sectors. Table 1 shows the total number of states and the 
number of pure states in each mass interval of fig. 1. We also show the expectation value 
of the number of partons averaged over all states in each mass interval. Evidently, a few 
low-lying states are pure, while there are almost no pure states among the high excitations. 
In fig. 2 we plot the spectrum of the odd glueballs for a; = and K — 2^, and in table 2 we 
quantify their purity. The lightest odd glueball contains, to a good approximation, only three 
partons, and the first excited state is almost purely a five-parton state. If we increase the 
mass of the quantum, the pair creation becomes somewhat suppressed. We have performed 
the diagonalization ior x — 1 and found, indeed, that the number of pure states increases. 
While the distributions of states for x = shown in figs. 1 and 2 are almost uniform within a 
band, for x = 1 there is an increase in the number of states with the average parton number 
near 2, 4, 6, 8, . . . for K = 24, and a similar increase in the number of states with the average 
parton number near 3, 5, 7, 9, . . . for K = 25. 

A striking property of figs. 1 and 2 is the rapid growth of the density of states with 
increasing mass. In fig. 3 we plot the logarithm of the number of odd glueballs vs. the mass 
for the data in table 2. For a certain range of masses the graph is approximately linear. 
The deviation from linearity for large enough mass is clearly due to the effects of the cut-off. 
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Our results indicate that the density of states grows roughly exponentially with the mass, 
exhibiting the Hagedorn behavior 

p(m) ~ m"e^"^ . (24) 

Thus, although the mass spectrum is discrete, it rapidly becomes virtually indistinguishable 
from a continuum. Prom our data we estimate that the inverse Hagedorn temperature is 
(3 (0.65 - 0.7) y/n/ig'^N). 

Another physical effect that is pronounced in our results (tables f and 2) is that the mass 
increases roughly linearly with the average number of partons. In fig. 4 we plot these results 
for x = and K — 2A (table 1). A heuristic explanation of this effect was suggested in [10] 
using the result of ref. [14] that the ground state energy of a system of n non-relativistic 
particles connected into a closed string by harmonic springs, which should be identified with 
M^, behaves as ~ for sufficiently large n. 

Since the low-lying states arc very pure, they can be well approximated by truncating 
the diagonalization to a single parton number sector. For instance, for ,t = and X = 24 the 
ground state has probability 0.999999 to consist of 2, 4 or 6 partons. Extrapolating to infinite 

we find 0.99996, and therefore this truncation is highly reliable. In table 3 we compare 
the full and truncated calculations. We have performed the truncated diagonalizations up 
to K — 40, and extrapolating these results to infinite K, we find the upper bounds f» 
A.SSg^N/n for x = and ISASg^N/n for x = 1. These are extremely close to the 

extrapolations from exact diagonalizations. This shows that, by judiciously truncating the 
space of states, certain eigenvalues can be determined to a good accuracy with relatively 
small diagonalizations. Similar approximations work well for the odd glucballs. For x = 
and K = 25 the lightest state has probability 0.999992 to consist of three or five partons. 
Extrapolating it to infinite K, we find 0.99994 . The advantage of this truncation is that 
we can access higher value of K (up to 51) than in the full diagonalization, and extrapolate 
more rehably. We find that the lowest eigenvalue extrapolated in this fashion is f» 
11.95 g'^N/n for x = 0, and 30.87 g'^N/n for x = 1. Good agreement of these 

values with extrapolations from exact diagonalizations gives us some confidence that our 
methods are consistent. We have used the Bulirsch-Stoer algorithm which has proved to 
be particularly efficient for extrapolating short series [15]. In fig. 5 we show the of the 
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ground state for x = as a function of 1/K. For < 24 we show the results of the full 
diagonalization, while for 26 < iiT < 40 we use the results of truncated diagonalization in 
the 2— ,4— and 6— bit sector. 

4. DISCUSSION 

Our calculations suggest that (l+l)-dimensional large- A?" QCD coupled to adjoint scalar 
matter captures some of the physical features of the higher-dimensional pure glue gauge 
theories. This model may also serve as a good test of the light-cone quantization methods 
that are promising to become a useful tool for studying the non-perturbative structure of 
the strong interactions. We mention below some interesting questions for further study. 

1. Even though the model is (1-1-1) -dimensional, we have observed a rich structure of 
glueball states. In fact, it is much richer than in the models coupled to matter in the 
fundamental representation of SU(A^), which have only one state per unit mass-squared [7, 
8]. It would be interesting to gain further insight into the Hagedorn behavior, eq. (24). 
Could there be exact formulae for the constants a and j31 Their determination is relevant 
to the properties of the deconfining phase transition. Such a transition has been recently 
studied in the model with adjoint fermions [11, 16]. 

2. We found that some low-lying states are exceedingly close to being eigenstates of the 
parton number. However, they fail to be exact eigenstates. We feel that this deserves a 
deeper understanding. If this property holds true in higher-dimensional gauge theories, it 
could lead to an enormously simplified picture of the low-lying states. 

3. The glueball states can be simply pictured as closed strings that cannot split or join 
in the IsiXge-N limit. Is there a continuum string description of our model? A study of 
this question seems to be a logical next step in the program of ref. [17], where a string 
representation of the pure glue theory was found. 
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Table 1. Numerical data for K = 24, x = shown in fig. 1; the 3.0 bin includes all states 
whose masses are 1.5 — 3.0; etc. 



M 


3.0 


4.5 


6.0 


7.5 


9.0 


10.5 


12.0 


13.5 


15.0 


16.5 


number of states 


1 


1 


1 


2 


7 


26 


103 


262 


598 


1013 


number of pure states 


1 


1 


1 


1 


2 


3 


4 


3 


1 


1 


a.A'orage ])ai't()ii iiuinl)ei' 


2.00 


1.00 


2.17 


1.98 


1.61 


1.96 


6.01 


7.36 


8.51 


9.91 



Table 2. Numerical data for K — 25, x — shown in fig. 2. 



M 


4.5 


6.0 


7.5 


9.0 


10.5 


12.0 


13.5 


15.0 


16.5 


18.0 


number of states 


1 


1 


1 


7 


26 


99 


277 


651 


1270 


1757 


number of pure states 


1 


1 


1 


3 


5 


5 


1 











average parton number 


3.00 


4.99 


3.09 


4.52 


5.12 


5.96 


7.18 


8.39 


9.69 


11.12 
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Table 3. Numerical values of the ground state mass for x = and x = 1 for even K. 





X = 


X = 1 


K 


M2 (full) 


M2 (2+4+6-bit) 


M2 (full) 


M2 (2+4+6-bit) 


12 


4.1064 


4.1064 


10.8712 


10.8712 


14 


4.1476 


4.1476 


11.0652 


11.0652 


16 


4.1789 


4.1789 


11.2219 


11.2219 


18 


4.2034 


4.2034 


11.3520 


11.3520 


20 


4.2230 


4.2230 


11.4624 


11.4624 


22 


1.2388 


1.2388 


11.5577 


11.5577 


24 


4.2518 


4.2519 


11.6411 


11.6411 












oo 


4.3 


4.33 


13.4 


13.43 



Table 4. Numerical values of the ground state mass for a; = and a; = 1 for odd K. 





X = 


X = 1 


K 


M2 (full) 


m2 (3+5-bit) 


M2 (full) 


m2 (3+5-bit) 


15 


10.5001 


10.5003 


26.0496 


26.0497 


17 


10.6569 


10.6572 


26.5323 


26.5326 


19 


10.7853 


10.7858 


26.9358 


26.9362 


21 


10.8925 


10.8932 


27.2797 


27.2802 


23 


10.9835 


10.9843 


27.5775 


27.5781 


25 


11.0616 


11.0626 


27.8387 


27.8393 












oo 


11.9 


11.95 


30.8 


30.87 
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FIGURE CAPTIONS 



Fig.l. The spectrum of even glueballs for K — 24, x — 0, M < 18; mass M is measured in 
units of ^/g^N/n and plotted vs. the expectation value of the parton number. 

Fig. 2. The spectrum of odd glueballs for K = 25, x = 0, M < 18. 

Fig.3. Logarithm of the density of states and a linear fit for K — 25 and x — 0. 

Fig. 4. Average number of partons as a function of mass for K = 24 and a; = 0. 

Fig. 5. Extrapolation towards infinite K of the ground state a,t x — 0; M^n/g^N is plotted 
versus 1/K. 
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